Generation of tunable Terahertz radiation using Josephson vortices 
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We propose how to control the THz radiation generated by fast moving Josephson vortices in 
spatially modulated (either along the c-axis or the o6-plane) samples of Bi2Sr 2 CaCu2 0g +( s and 
related superconducting compounds. We show that the JVs moving in a subset of weaker junctions 
can generate out- of-ab -plane and outside-the-cone Cherenkov radiation. The ab-plane modulation 
of superconducting properties (achieved, for instance, by ion irradiation lithography) can result 
in transition radiation within certain frequency windows, i.e., allowing the design of tunable THz 
emitters and THz photonic crystals. 

PACS numbers: 74.25.Qt, 74.72.Hs, 41.60.-m 



The recent growing interest in terahertz (THz) science 
and technology is due to its many important applications 
in physics, astronomy, chemistry, biology, and medicine, 
including THz imaging, spectroscopy, tomography, med- 
ical diagnosis, health monitoring, environmental control, 
as well as chemical and biological identification Q. The 
THz gap, that is still hardly reachable for both elec- 
tronic and optical devices, covers temperatures of biolog- 
ical processes and a substantial fraction of the luminosity 
remanent from the Big Bang [lj. 

High-temperature Bi 2 Sr2CaCu20 8 +5 superconductors 
have a layered structure that allows the propagation of 
electromagnetic waves (called Josephson plasma oscilla- 
tions 0, EJDI 0- ) with Josephson plasma frequency ui j . 
This is drastically different from the strong damping of 
electromagnetic waves in low temperature superconduc- 
tors. It has been recognized (see, e.g., that the 
Josephson plasma frequency lies in THz range. A possi- 
ble way to generate THz radiation in Bi2Sr2CaCu20s+5 
and related compounds is to apply an in-plane magnetic 
field H ao an d an external current J|| c perpendicular to the 
superconducting layers (i.e., along the c axis). Josephson 
vortices (JVs) induced by H ao and driven fast by the o 
axis current emit THz radiation (e.g., jjj,|9j). However, 
it was shown 0, 0, 0] that the radiation propagates 
only along the plane of motion of the JVs and decays 
in the c-direction. This strong confinement of THz radi- 
ation, its rather restricted controllability, and the large 
required c-axis current, all limit potential applications. 

Recent developments in sample fabrication can pro- 
duce superconductors with alternating layer s 1131 , as well 
as spatially-modulated both pinning (e.g., [14() and su- 
perconducting coupling between layers. We propose em- 
ploying these capabilities to create modulated structures 
for a new generation of THz emitters based on moving 
JVs. 



Radiation from a weak junction: A subset of weaker in- 
trinsic Josephson junctions in Bi2Sr2CaCu20s+5- based 
samples can be made using either (i) the controllable in- 
tercalation technique ILJ, (ii) Chemical Vapor Deposi- 
tion (CVD) (see, e.g., (lfj), or (iii) via the admixture of 
Bi 2 Sr 2 Cu20 6+< 5 and Bi 2 Sr 2 Ca2Cu30io+5 0- The moti- 
vation for making such artificial samples (Fig. lb) is that 
the radiation produced by the usual Bi2Sr2CaCu20s+i 
materials is confined to the ab-plane. In other words, if 
all planes (i.e., all junctions) are identical, the Cherenkov 

radiation generated by fast moving jvs 0, EE El E3 

propagates only along the direction of the JV motion. 
Electromagnetic (EM) waves along the c-direction can- 
not be generated in identical-planes samples because the 
maximum velocity of the JVs is smaller than, or of the 
order of, the smallest velocity of the c-axis EM waves. 
Thus, previous proposals for pseudo-Cherenkov radia- 
tion 0, EEl G2 in standard Bi 2 Sr 2 CaCu20 8+5 sam- 
ples do not refer to the usual Cherenkov radiation of a 
fast relativistic particle because the radiation generated 
by those JVs is narrowly-confined, i.e., it does not prop- 
agate within a Cherenkov cone. We propose using an 
artificial stack of thin Bi2Sr2CaCu20s+5 sheets having 
several Josephson junctions which are weaker than the 
other intrinsic ones. The JVs in the weaker junctions can 
move much faster than the Josephson plasma waves in 
Bi2Sr2CaCu208+,5 materials, producing electromagnetic 
waves propagating both parallel and perpendicular to the 
CuOi layers (Fig. la,b). Changing the c-axis current 
changes both the JV speed and the radiation power. 

Transition radiation: A completely different design 
involves samples with periodically modulated 01 ( m 
the ab-plane) superconducting properties which are now 
uniform along the c-axis. For example, spatial (in- 
plane) variations of the Josephson maximum c-axis cur- 
rent J c can be obtained by using irradiation of a standard 
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FIG. 1: Cherenkov radiation generated by a fast Joseph- 
son vortex (located at x = Vt) moving in a weaker junc- 
tion, (a) Magnetic field distribution H(x — Vt,y) in units of 
$o/2iv\ c \ ab for J?/J c = 0.2, s w e/(se w ) = 1.2, V/V™^ = 0.9. 
The "r unning" c oordinate, x — Vt, is measured in units of 
7s/(7Ty / v 2 /3 2 — 1), while the out-of-plane coordinate y is nor- 
malized by s/(ttvV/3 2 - 1), where /3 = tt J c es w /2 J™e w s. The 
moving vortex emits radiation propagating forward. This ra- 
diation forms a cone determined by the vortex velocity V. 
(b) A suggested experimental set up: in a weaker junction 
an out-of-plane current J\\ c drives a Josephson vortex with 
velocity V, which is higher than the minimum velocity c m i n 
of the propagating electromagnetic waves. 



Bi2Sr2CaCu20g+5 sample covered by a modulated mask 
(see, e.g., 14]). For such a sample, the Josephson vor- 
tices can produce transition Terahertz radiation (Fig. 2a) 
even if the vortex velocity is slower than c m ; n . The termi- 
nology "transition radiation" 0] refers to the radiation 
produced when the vortex alternatively transits from one 
"layered medium" to another. The spatial periodicity in 
the a6-plane controls the radiation frequency. 

Josephson vortex as a fast relativistic particle. — Be- 
sides potential applications to THz technology, the mo- 
tion of a JV in layered superconducting materials has 
significant scientific interest because it can mimic some 
properties of fast relativistic particles. Indeed, it is well- 
known that the Sine-Gordon equation, describing the mo- 
tion of a JV in a conventional Josephson junction, is in- 
variant under a Lorentz transform, where the speed of 
light is replaced by the Swihart velocity [2(J . For a stan- 



dard Josephson junction, the Swihart velocity restricts 
both the maximum speed of small magnetic field per- 
turbations and the maximum vortex velocity. In layered 
structures, the equation describing the Josephson vortex 
and its magnetic field H{x,y,t) becomes nonlocal |2l| 
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Here, the Josephson plasma frequency ojj is ~ 1 THz, 
depending on doping and temperature; the two length 
scales A c = j\ ab and A Qb = 2000/^/1 - T 2 /T 2 A, with 
T c = 90K and 7 ~ 300, determine the characteristic 
scales of magnetic field variations parallel (along the x- 
axis) and perpendicular (along the y-axis) to the super- 
conducting layers (Fig. lb). Eq. describes the prop- 
agation of the gauge invariant phase 4> = Xi ~ X2 + 
2irA y s/&Q along the junction where the Josephson vor- 
tex moves, where \i an d Xi are the phases of the wave 
functions of the superconducting condensate of the Cu02 
layers forming the junction, and A y is the y-component of 
the vector potential. Here, s denotes the thickness of this 
junction and $0 is the magnetic flux quantum. Eq. (J2J 
has its own space and time scales: I and u>~ 1 . There Kq 
denotes the modified Bessel function usually employed for 
the magnetic field distribution in superconductors. When 
all Josephson junctions are the same (e.g., as in standard 
Bi2Sr2CaCu208+5 samples), uj c = wj , s — s, and the size 
I of the soliton at V — is defined as a product of the 
anisotropy parameter 7 and the spacing s between CuC>2 
planes (I ~ 7s). For a weaker Josephson junction, the 
Josephson soliton is more elongated, I — l w — "/sJ c /J™, 
and 

(jJ c _ u w is 

determined by the parameters of the 
weaker junction: ui w = u>jy/ ' s w J™e/ sJ c e w with J c , e 
are the maximum superconducting current and dielec- 
tric constant of the intrinsic junctions; the index w refers 
to the weaker Josephson junction. 

Eq. Q for the magnetic field H allows us to obtain 
the spectrum of EM waves propagating in the layered 
structure 
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of the propagating EM wave are determined by three 
standard conditions: (i) Eq. Q, (ii) uj = k x V, and (iii) 
the minimum wavelength k y < n/s. However, the JVs 
cannot move above the maximum speed V„ lax determined 
by the nonlocal Eq. J2J): V max ~ ujj^s/2 for identical- 
junctions sample and V|™ ax ss lu w I w for samples with some 
weaker junctions. Therefore, the radiation propagating 
within the Chcrenkov cone with 9 =/= cannot be gen- 
erated using standard Bi2Sr2CaCu20s+5 samples, while 
it can be achieved using the proposed artificial stack 
or c-axis-modulated Bi2Sr2CaCu20s+5 samples. Using 
Eqs. we calculate perturbatively the Cherenkov ra- 

diation (see Fig. la): 
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where the dynamical-soliton size l(V) decays when 
the vortex velocity V 0] increases as: l(V) = 

i (i + Vi-v-7(K£j 2 )/2- 



l a [1 + sjl - V 2 /(V™ ax ) 2 )/2. The y-component k y of 
the wave vector is determined by the dispersion rela- 

Due to the 
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rather unusual dispersion relation i.e., the decrease of 
k y {k Xl u) = k x V) with increasing k x , the electromagnetic 
waves are located outside the Cherenkov cone (Fig. la), 
which is drastically different from the Cherenkov radia- 
tion of a fast relativistic particle. The new type of ra- 
diation predicted here could be called outside-the-cone 
Cherenkov radiation. 

Zone structure of the transition radiation. — Here 
we consider standard Bi2Sr2CaCu20s+<5 materials with 
identical junctions. As mentioned above, previousl y p ro- 
posed pseudo-Cherenkov-like radiation pj UlA liil LL3 is 
confined to the ab-plane, and requires a rather strong 
c-axis current. In contrast to this, we consider an alter- 
native approach to generate THz radiation using stan- 
dard Bi2Sr2CaCu208+5 materials which allows to control 
the spectrum of the radiation. This proposal might be 
the easiest to experimentally implement, among the ones 
listed here, since it does not require samples with weaker 
junctions. 

The modulation of superconducting properties along 
the ab-plane could be done by using, e.g., either lithog- 
raphy or by controlling the JV dynamics via their interac- 
tions with other type of vortices (pancake vortices; with 
its density determined by the out-of-plane field H c ). In- 
deed, the c-axis field, H c , creates another type of vortices, 
called pancake vortex stacks, which have slow dynamics 
and attract Josephson vortices |U • The interaction 
of these two types of vortices changes the dynamics of the 
JVs 0, Hj| and affects their motion in a similar manner 
as spatially modulating the critical-current with a period 
which can be easily tuned by the out-of-plane field H c . 
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FIG. 2: Transition radiation emitted by a Josephson vor- 
tex (located at x = 0) moving through a spatially modu- 
lated (along the ab-plane) layered superconducting sample, 
(a) Magnetic field distribution H(x,y) (at a certain time mo- 
ment, say t = 0) in units of fj,i^os/2Tv\^ b a vllT for V/c m i n = 0.8, 
a val /l(y = 0) = 1. The in-plane and out-of-plane coordinates 
x and y are normalized by the core size I of a static Joseph- 
son vortex and 2s, respectively. (b,c) The i-component k x of 
the wave-vector of the radiation versus frequency. In contrast 
to the Cherenkov radiation (Fig. 1), the phase velocity k x /uJ 
of the transition radiation could be positive or negative (b), 
resulting in waves propagating both forward and backward 
with respect to the Josephson vortex motion. The radiation 
frequency has forbidden zones, shown by red strips in (c), 
when the vortex moves relatively slow. This suggests the re- 
markable possibility of THz photonic crystals in modulated 
layered superconductors, which is potentially important for 
applications. 



To simplify the problem, we study the radiation of 
Josephson vortices with modulated u c (x). Note that the 
modulation of any other parameters (e.g., l(x) or both 
l(x) and u) c (x)) in Eq. @ results in the same qualitative 
results. For periodically-varying properties (with spatial 
period a va r) of intrinsic junctions, the frequency w and 
the in-plane wave vector k x of the emitted EM waves 
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are related by the condition lj — (k x + 2irm/a vaT )V , 
with m = 1,2,.... This is the momentum conser- 
vation law: the momentum of a moving JV can be 
transferred not only to the EM waves but also to the 
modulated medium, in full analogy with electron mo- 
tion in a periodic potential. Following this analogy, 
we derived the "zone structure" of the radiated waves 
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Here k x (m) is related to the momentum 27rm/a var trans- 
ferred to the modulated medium by k x = Lu/V—2irm/a var 
and /i m is the Fourier component of 2[1 — uj c (x)/ujj]. For 
a chosen zone-number m, the radiation is emitted in a 
certain frequency window w m i n (m) < uj < w max (m), with 



y mm, max 



(m) 



2mn c n 



V 2 



V 2 



V 



J \ 2nmc n 



(c min ) 2 - V 2 
V 2 



(8) 



Since the frequencies w m i nj max change with changing a var , 
the generation of THz radiation in a certain frequency 
range could be designed by adjusting the period of the 
spatial modulations of ion irradiation during the sample 
preparation. When the out-of-plane magnetic field H c is 
applied, the transition THz radiation might be tuned by 
varying H c . 

In contrast to the outside-the-cone Cherenkov-like 
radiation (Fig. la), our proposed transition radiation 
(Fig. 2a) propagates both forward and backward in space. 
Indeed, the in-plane component k x of the wave vector 
and, thus, the corresponding phase velocity both change 
sign (Fig. 2b, c). Also, the EM waves running backward 
can be directly seen from the magnetic field distribution 
shown in Fig. 2a. For fast Josephson vortices moving 
with velocity close to c mm , the frequency zones over- 
lap for different zone number m (Fig. 2b). However, for 
slower speeds, we obtain the forbidden frequency ranges 
w max (m) < uj < u> m - ln (m + 1) (see, Fig. 2c) of radi- 
ated electromagnetic waves. Thus, such materials might 
be used for making tunable THz photonic crystals con- 
trolled, e.g., by the changing out-of-plane magnetic field 

Conclusions. — A grand challenge is to controllably 
generate electromagnetic waves in Bi 2 Sr 2 CaCu208+5 and 
other layered superconducting compounds because of its 
Terahertz frequency range [3, la M ■ Considering recent 
advances in sample fabrication |l3l Hif , we propose sev- 
eral experimentally realizable systems for generating con- 
tinuous radiation in a controllable frequency range. 
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